Scenario for Fractional Quantum Hall Effect in Bulk Isotropic Materials 
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We investigate the possibility of a strongly correlated Fractional Quantum Hall (FQH) state in 
bulk three dimensional isotropic (not layered) materials. We find that a FQH state can exist at 
low densities only if it is accompanied by a staging transition in which the electrons re-organize 
themselves in layers, perpendicular to the magnetic field, at distances of order the magnetic length 
apart. The Hartree energy associated to the staging transition is off-set by the correlation Fock 
energy of the 3D FQH state. We obtain the phase diagram of bulk electrons in a magnetic field 
subject to Coulomb interactions as a function of carrier density and lattice constant. At very low 
densities, the 3D FQH state exhibits a transition to a 3D Wigner crystal state stabilized by phonon 
correlations. 

PACS numbers: 73.43.f, 11.25.Hf 
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The Quantum Hall effect is intimately linked to the 
low-dimensionality of the sample and the fact that the 
magnetic field quenches the kinetic energy of the electron 
liquid in 2 dimensions (2D). The Quantum Hall effect has 
also been observed in 3 dimensions (3D) in the Bechgaard 
salts [l|, Q but only at integer filling factors and only 
in layered samples. The 3D Quantum Hall states can 
be explained as a series of 2D states, weakly coupled so 
that the band-width in the z-direction be smaller than 
the 2D Landau gap. So far, 3D FQH states have not 
been observed, even in layered samples, due to their low 
electron mobility and small many-body gap. 

Recent experimental interest in graphene, graphite and 
Bismuth has focused on the Dirac nature of carriers and 
on the Integer Quantum Hall (IQH) effect. The Dirac 
dispersion in these materials gives electrons large mean- 
free paths and large Landau gaps. As a result, large IQH 
plateaus are observed, even at room temperature. A 3D 
IQH effect has been predicted in graphite 3] ■ Band struc- 
ture considerations suggest that, at sufficiently strong 
magnetic field, the n = and n — 1 Landau levels of 
graphite are gapped and a single 3D IQHE plateau ex- 
ists. The layered structure of graphite is essential to the 
existence of the 3D IQH. Recent experiments in Bismuth 
suggest the possibility of a 3D bulk strongly correlated 
electron state in the fractional filling regime ij] ■ The ex- 
periment reports quasi-plateaus in p^y at fields which are 
integer times larger than the field required to drive the 
system in the quantum limit. 

These results are unexpected since Bismuth, unlike 
Graphite, is an isotropic material. There are only two 
known scenarios for the existence of a Quantum Hall ef- 
fect in 3D: it can be either a band-structure effect, such 
as in graphite, or an interaction effect, where a 2kp in- 
stability opens a gap at the one dimensional Fermi level 
in the direction of the magnetic field. The Fermi level 
is then pinned in the many-body gap and the system 
exhibits a Quantum Hall effect. However, by construc- 



tion, both these scenarios lead to an integer, and not 
fractional state. While it is easy to generalize the Laugh- 
lin [5] and Halperin ^ states to a 3D FQH state 0,11 
with correlations between layers, they were never found 
to be the variational ground-states in an isotropic ma- 
terial. Contrarily, it was found that the magneto- 
plasmon gap closes and a phase-transition to a crystal 
state occurs whenever the distance between the layers is 
roughly smaller than half the magnetic length. 

Previous work on quantum well bi-layer and multi- 
layer systems suggests several possible electronic arrange- 
ments for a 3D material in a magnetic field: MacDon- 
ald and coauthors considered the energetics of 3D 
Halperin states, which are adaptations of the Laughlin 
wave functions to ensure that electrons in adjacent lay- 
ers maximize their separations. At very small inter-layer 
separations, they found these states crystalize due to a 
collapse of the magnetoroton gap. MacDonald [3] also 
established that in multilayers, in the Hartree Fock ap- 
proximation, the gain in exchange energy from distribut- 
ing electrons unequally between the layers can exceed 
the electrostatic cost of increased inter-layer Coulomb 
energy leading to a staging transition. Ref. [lo| pro- 
posed a third type of candidate state, the spontaneous 
inter-layer coherent miniband state, consisting of a com- 
bination of lowest Landau level (LLL) states in each layer 
which forms a band in kz- However, no state exhibiting 
a 3D FQH effect has been found as the ground-state of 
an isotropic material in B field. 

In this paper, we revisit these issues and try to address 
the question of whether a FQH effect could exist in the 
3D compound. We propose a new variational ground- 
state, a correlated 3D FQH-staged state in an isotropic 
material. We compare four different potential ground 
states for materials in B field in the ultra-quantum limit 
with strong Coulomb interactions and find that 3D FQH 
ground-states do not occur in isotropic materials unless 
they are accompanied, at very low densities, by stag- 
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ing transitions similar to the ones observed in graphite 
intercalation compounds [ll| and predicted to occur in 
multi-quantum- well systems by MacDonald In the 
parameter range relevant to both graphite and Bismuth, 
we obtain the phase diagram as a function of the electron 
filling fraction, the ratio of inter-layer separation to mag- 
netic length, and the magnitude of the c-axis hopping. At 
high carrier density and zero c-axis hopping, a (staged) 
integer quantum Hall liquid is energetically favorable. As 
the filling in occupied layers decreases below ly = 1, liquid 
states with correlations between layers can develop, and 
become energetically favored because of their decreased 
inter-planar Coulomb energy. At extremely low densities, 
staged crystal states tend to be energetically favored. A 
SILC miniband state becomes favored over the staged in- 
teger liquids, but not their fractional counterparts, at the 
c-axis hopping relevant to graphite. Thus we find a pa- 
rameter regime in which fractionally filled, staged quan- 
tum Hall liquid states are the expected ground states of 
the layered system. 

Let us consider an isotropic material subject to a mag- 
netic field parallel to one of its crystallographic axes. 
The magnetic field quenches the kinetic energy of elec- 
trons in lattice planes (layers) perpendicular to the mag- 
netic field. A plausible scenario, borrowed from MacDon- 
ald 's work on staged states in multiple quantum wells 
is that the relatively small inter-layer separations in 
isotropic materials favor staged states which consist of n 
layers of electron gas of density — Sa, and one layer 
of density We + nda, where CTe is the initial electron den- 
sity per layer, and Sa is the staging density. For crys- 
talline states, staging increases the distance between elec- 
trons which lie directly above one another, decreasing 
the crystal's Coulomb energy at small inter-layer sepa- 
rations; for liquid states staging becomes advantageous 
when the gain in Fock energy by increasing the electron 
density in the filled layers outweighs the Hartree cost of 
distributing charge unevenly between the layers. Thus 
as the separation between layers decreases, staged states 
become increasingly energetically favorable. 

We calculate the staging Hartree energy as described 
in 0]. We model the system as a stack of planes, each of 
uniform charge density cr^, where ai — —nda if i|(n -I- 1) 
and 6a otherwise. This charge density is comprised of 
an immobile positive background (ionic) charge density 
ae in each plane, and a staged electronic charge density 
—ae + ai. n is called the staging number. The 3-d charge 
density is conveniently described in Fourier space as a 
sum of Kroeneker 6 terms 
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where V is the total sample volume, and d is the lat- 
tice constant in the direction of the magnetic field. This 
charge density configuration has a Hartree energy Eg — 
U|p(g,)|2,or: 



eV 



Es = 



{5vY d 



el 



121 



{n^ + 2n) 



(2) 



where vq = 2TTl^ae, and = 2TTp6a is the electronic 
filling staged from each layer. One could also consider 
multiply staged states, in which three or more different 
staging numbers n exist. However, we numerically found 
that the states of lowest energy are those with charge 
distributions of the form in Eqn. ([T]), described by a 
single staging number. 

The first class of candidate states that we consider are 
the staged liquid states. The simplest such staged liquids 
consist of n layers depleted of electrons, and one layer of 
electronic charge density ae + nSa (which we will call the 
occupied layer) in a quantum Hall liquid state. The total 
energy of a staged liquid state is then: 

n 1 
E=—-Ei{vo-Sv) + ^—Ei{iyo + nSiy) + E, (3) 
n + 1 n + I 

where n is the staging number, Ei{u) is the energy of 
the liquid at filling u, Eg is the Hartree energy given by 
([2]), z^o is the mean filling factor, and Q < 5v < v^^ is the 
amount of charge staged out of each layer. 

As Laughlin states have good in-plane correlation ener- 
gies, we consider first the case where the occupied layers 
are in a Laughlin state. We will refer to Laughlin states 
at filling fraction 1/m in the occupied layers as (0,m,0) 
states. Since the filling fraction in these layers is then 
fixed at a value oi v = l/m for some odd integer to, the 
quantities v, vq and n are related by n + 1 = ^ . Using 
the optimal staged liquid state for a given vq and d/l 
can be found by optimizing over Sv and n. We find that 
optimal states in the regimes of interest are fully staged 
{Sv = I'o), so that the first term of the right-hand side of 
([3|) vanishes. 

At small d/l, the simple uncorrelated liquid state de- 
scribed above gives way to liquids with inter-layer corre- 
lations |9|, which we call Halperin liquid states. These 
states have wave functions of the form 
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where k indexes the layer, and i,j index particles within a 
layer. We refer to states of the form @ as (toi,TO2,toi) 
states. Fig. [T] shows the energies for several of these 
states as a function of inter-layer separation and density, 
calculated using Monte Carlo techniques (see Appendix) . 

When d/l is very small, one can imagine a low en- 
ergy ground state of staged Halperin liquid states, with 
n layers completely depleted of electronic charge for ev- 
ery occupied layer at filling v = 1/(2toi -I- TO2). In the 
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layer repulsions favor an in-plane lattice that is square as 
this permits larger separations between sites in neighbor- 
ing planes As d decreases, multiple phase transitions 
between different stackings of squares in the vertical di- 
rection occur; reference ^] identifies states in which elec- 
trons at the same in-plane co-ordinates are separated by 
M = 2, 4 or 5 layers, as well as states with incommensu- 
rate stackings. In this work we consider only true crystal 
states with well-defined spatial periodicity, and will not 
include structures with incommensurate stackings. 

Numerical studies 1^ indicate that the total energy 
of these crystalline states is well approximated by 



FIG. 1: Energies of unstaged Laughlin (dashed lines) and 
Halperin (solid lines) liquid states as a function of inter-layer 
separation over magnetic length (d/l). All energies are re- 
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ported as energies per electron, in units of where I is 
the magnetic length. At large separation the energy increases 
with d/l as the correlations in Eqn. Q do not optimally min- 
imize the intra-planar Coulomb interactions and their energy 
grows (see text). The minimum energy occurs when corre- 
lations effectively separate electrons from both their in-plane 
and out-of-plane neighbors. 



layered system with a given mean filling per layer vq , the 
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staging number is fixed by n 4- 
in turn fixes the separation (n 



i/o(2mi+m2)' 

l)d between occupied 
layers, and the energy of such a state is 



Ei{{n + l)d) + Esid,n,v) 
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where Ei{{n + l)d) is the Coulomb energy of the state 
(j4|) at inter-layer separation (n -1- l)d, extrapolated from 
the Monte Carlo results of fig. [U and Eg is the Hartree 
staging term. Fig. [1] shows an important property of 
the energies as a function of d/l. The wave function (|4]) 
contains correlations which separate electrons in adja- 
cent layers, at the expense of decreasing their average 
separation in a given layer. Since these correlations do 
not depend on d/l (the out of plane correlations in Eqn. 
(|^ depend on z^'^'' — z^'^ ^\ i.e. only on the in-plane 
distance between particles in different layers), as d/l in- 
creases the energy of the Halperin states also tends to 
increase. The inter-planar Coulomb interaction becomes 
weaker, so that the energetic gain from inter-layer cor- 
relations decreases, while the energetic cost of decreased 
separation in-plane relative to the Laughlin state of the 
same filling remains constant. At large d/l the Laughlin 
liquid becomes the ground state. 

In two dimensions, it is well known that for a low den- 
sity of carriers, a transition occurs from a FQH state 
to a triangular lattice Wigner crystal of lattice constant 

— I I vAE j xhe Wigner crystal is energetically fa- 
vored over the FQH liquid state at filhngs ly < 1/7 [l^ . 
In multi-layer systems at small d {d/l < 0.9 J inter- 
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Here Eq is the classical contribution, given by the energy 
of the Madelung sum, E2 is the phonon contribution, and 
£"3 is an extra higher moment contribution calculated by 
a fit to a numerical evaluation of the Coulomb energy. For 
the low fillings at which Wigner crystals are expected to 
occur, the total energy is well approximated by the first 
two terms. 

The classical contribution Eq can be calculated us- 
ing the Ewald method, described in detail in 13]. We 
slightly modify this method to account for the fact that 
the background charge is localized in the planes. At each 
value of d, energies were tabulated for a variety of crys- 
tal stackings and stagings; the configurations of lowest 
energy were selected. We estimate Ej, by numerically 
extrapolating the results of 9] as a function of nd/l to 
the values of interest here ((n -I- l)d/l « 0.2 for most 
optimally staged crystal states). 

We compute the contribution of lattice vibrations E2 
by postulating a Lam-Girvin [l^ variational wave func- 
tion of the form 



exp n 



(7) 



Here Rj = {Xj,Yj) is the position of the lattice site, 
= [xj — Xj) -\- i{yj — Yj) is the deviation of the parti- 
cle's position from this site in two-dimensional complex- 
ified coordinates, and <pji(f) the lowest Landau level 
coherent states, 
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The correlation matrix B is a variational parameter; min- 
imizing the variational energy gives 
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FIG. 2: Calculated energies (per electron) of the unstaged 
crystalline states, as a function of d/a where a is the crystal 
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lattice constant. Energies are shown in units of Top: 
Coulomb sum contribution (i?o). As d/a decreases several 
phase transitions occur to crystal stacking patterns with a 
larger unit cell in the z direction, consistent with the findings 
of Bottom: Phonon contribution (£^2) for the crystal 
states of lowest Eq. As d/a decreases the mean inter-electron 
distance decreases and E2 grows. 



FIG. 3: A comparison of the energies of the (0,1,0) and 

(1, 1, 1) staged liquids with those of SILC miniband states 

as a function of mean filling for d/l = 0.1. Energies and 
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t± are shown in units of At the physical value t± = 
0.033-e^/eZ the SILC miniband state is energetically favored at 
high densities, and has lower energy than the staged integer- 
filled liquid state. The vq dependence of the liquid energies 
can be understood from fig.[l]by noting that at fixed d/l and 
occupied filling u, nd/l — vd/luo). 



where ujl and lot are the longitudinal and transverse 
phonon frequencies, respectively, and E is the energy of 
the phonon modes. The phonon frequencies are calcu- 
lated using the method of For multiple-site unit 
cells, we fix as a single function independent of site 
indices within the unit cell. In this case the energy is 
given by ([9]), with lut,ujl replaced by the square roots 
of the eigenvalues of the 2x2 matrix formed by averag- 
ing the second order correction to the Coulomb potential 
over all sites in the unit cell. 

Fig. [5] shows the energies Eq and E2 of the crystalline 
states over a range of inter-layer separations. At very 
small d/a, the dominant interaction for a M-layer crys- 
tal stacking is between an electron and its translates M 
layers above and below; the associated energy increases 
roughly as -pg. To avoid this cost, the crystal under- 
goes several transitions to stackings of higher M over 
the range of inter-layer separations shown in fig. [5] @; 
at sufficiently small d/a even the M = 5 stacking does 
not adequately separate charges from their closest verti- 
cal neighbors and the energy grows rapidly. Though we 
computed the energy of all three crystal structures (2- 4- 
and 5- layer crystal stackings) at each value of d/a, only 
the lowest of these energies is shown in the fig. [2] 

At the very small values oi d/l found in graphite and 
bismuth, none of these crystal structures represent a sta- 
ble ground state. We therefore allow for Wigner crystals 
with staged charge densities, requiring that charge de- 
pleted layers be completely emptied of electrons {5a = 
CTe)- The total energy for the staged Wigner crystal is 
given by 

E = E,,y, + Es (10) 



where i?crys is given by ([6]), using the values of Eq and 
E2 shown in fig. [2] at the effective inter-layer separation 
(n -f l)d, and Eg is given by ([2]) with 5v = lyQ. At a 
given inter- layer separation d and mean filling i^q, we find 
the energetically optimal Wigner crystal by choosing the 
staging number which minimizes (|10p . 

The final candidate state we consider is the Sponta- 
neous Inter-Layer Coherent State (SILC) miniband of 
Hanna, Dfaz-Velez, and MacDonald 10]. The SILC is 
given by the Slater determinant state 

W -11^1x10) (11) 
where the creation operator generates the state 

(r1C],^|0) = J2e^^''%,x{r^ (12) 

V^^p j=i 

^,x(r) = -A_^(,„,rf)e^x./'%-(-^)V2.^ 

Here, X is the center of the Landau strip, and x{z — jd) 
describes the vertical localization of the electron to the 
j"^ plane; Np is the number of planes (layers). The prod- 
uct X is over all states in the lowest Landau level (LLL), 
but only the lowest energy states in the momentum band 
are filled: —nv/d < q < TTv/d. 

The energy of the SILC is given by [§| 

Ec i^e^ f d^r -r"^ jii^ /sin(7rjV)\^ 2tj^ sin(7rz/) 

TVe ~ el ^ J 47r/2 Y TTjiy J TTiy 

(13) 
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FIG. 4: Approximate phase portrait for graphite in the range 
B ~ SOT. SILC miniband states are navy blue; the Laughlin 
liquids are colored dark blue ((0,3,0) state), pale blue ((0,5,0) 
state), and turquoise ((0,7,0) state). The Halperin liquids are 
colored orange ((1,1,1) state), red ((1,3,1) state), and brown 
((1,5,1) state). Crystalline states are shown in green. 



where t±_ is the hopping matrix element in the stack- 
ing direction. Fig. [3] compares the energies of the SILC 
miniband states at d = 0.1 for several t± values with 
the energies of the staged liquid states. The energy of 
the SILC miniband state depends strongly on the c-axis 
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hopping matrix element t±. For t± = 0.01||-, at the val- 
ues of d pertinent to graphite the SILC miniband state 
is never a ground state, as shown in the fig. [3] At in- 

2 

termediate values { ^ 0.03||-), the SILC miniband 
state is the ground state for sufficiently high densities. 

2 

For > 0.05||-, the SILC miniband state is the ground 
state for all densities shown. 

We now compare the energies of the three types of 
states outlined before in a parameter regime experimen- 
tally attainable for graphite. We model graphite as a 
quantum Hall multilayer of graphene sheets. Within each 
graphene plane, wc assume that the magnetic length is 
sufficiently large that the positive charge of the graphene 
crystal is well approximated as a uniform surface charge 
We- Further, we take the inter- plane separation d to be 
fixed, with a value of approximately d — 3.4 A. We con- 
sider relatively strong magnetic fields {B « 30 T). 

In graphite, at i? = 30 T, the relevant length scales 
IS I = 47 A, giving d/l « 0.07. The c-axis hopping in 
the LLL, as calculated from the c-axis bandwidth of the 
LLL 01, is i_L ~ lOmeV and the Landau gap is given by 



30 meV. The energy scale is set by 



0.31 eV. 



We assume that the mean filling fraction can be tuned 
independently of the magnetic field. 

The expected phase portrait for the system is shown 
in fig. m At constant d/l, as the mean filling decreases 
(in real materials, these would be differently doped sam- 
ples) , the ground state shifts from a SILC miniband state 
through (1,1,1), (0,3,0), (1,3,1), (0,5,0), (1,5,1), and 
(0, 7, 0) liquid states; at sufficiently low vq a phase transi- 



tion to a staged crystalline state (shown in green) occurs. 
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At physically relevant values of t± « 0.03^, the SILC 
miniband state is the ground state only at high densi- 
ties where integer-filled Laughlin liquids would otherwise 
occur. 

The structure of the phase portrait reflects the com- 
petition between Hartree and Fock energy contributions. 
The Fock energy decreases with increasing density, and 
thus always favors maximal staging. As vq decreases at 
fixed d, the Fock energy of staging to a given filling u 
remains fixed, while the Hartree energy due to staging 
increases approximately linearly in n + 1 — v /vq. Hence 
a series of transitions to states of lower filling in the occu- 
pied layers occurs. Once v < 1/7, the crystalline states 
have lowest energy. At fixed a transition occurs be- 
tween a Halperin state at higher (and hence smaller 
separation nd between occupied layers) and a Laughlin 
liquid at lower i/q (favored at larger nd). The effect of 
increasing d at fixed is similar: it increases the effec- 
tive cost of staging, and hence favors less filled states; at 
fixed V a phase transition between Halperin and Laughlin 
liquid states occurs. 

The crystalline states investigated in this paper have 
all assumed 'pancake' charge distributions confined to 
individual lattice planes. However, c-axis hopping will 
allow the electron charge distribution to spread out in 
the z-direction. For weak hopping, the crystalline states 
should be stable throughout much of their phase dia- 
gram, but eventually as increases, states like the SILC 
are preferred. In the continuum, one can define a fam- 
ily of variational three-dimensional crystalline states as 
generalizations of the two-dimensional quantum Wigner 
crystal of Maki and Zotos [l3| , writing 



* = det 



(14) 



where is the electron position projected onto the (x, y) 
plane, and z is the c-axis coordinate, and where 



(15) 



is a product of the in-plane lowest Landau level coher- 
ent state and a trial single particle wavefunction describ- 
ing the localization of the electron along the c-axis. An 
obvious choice for ip{z) would be a harmonic oscillator 
wavefunction with a length scale A that is a variational 
parameter, determined by the local curvature of the self- 
consistent crystalline potential. In this paper, we under- 
estimate the stability of our crystalline states; allowing 
the electrons to delocalize somewhat in the transverse 
dimension will lower their energy. We will report on cal- 
culations based on these states, and their obvious ex- 
tensions to correlated Wigner crystal states described in 
Eqn. [71 in a future publication. 

We conclude that fractionally filled quantum Hall 
states cannot occur as stable ground-states in isotropic 
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materials under realistic B fields for which the magnetic 
length is much larger than the c-axis lattice constant un- 
less the density is sufficiently low so that a staged FQH 
state becomes energetically favorable. At even smaller 
densities, the FQH staged ground-states give way to a 
staged 3D Wigner crystal with large unit cells in the di- 
rection parallel to the magnetic field. 
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Appendix: To use the Monte Carlo method initially used 
by Laughlin for FQH systems [11 in a layered system, 
we chose periodic boundary conditions and compute the 
Coulomb energy of each configuration as a Ewald sum 
over repeated copies of the fundamental cell. Particle 
configurations are generated with a probability dictated 
by the squared wave function; the energy is calculated 
by averaging the Coulomb potential over a large number 
of configurations. The fundamental cell is a square in 
plane of side length ^J^h^WN^s^ where iV2D is the to- 
tal number of particles per layer. These dimensions are 
chosen to enclose the maximum possible area over which 
the distribution of electrons in a droplet is essentially 
uniform, minimizing boundary effects. At each step in 
the simulation, only particles which are found within the 
fundamental cell are included in calculating the energy, 
though particles outside this region may return at a later 
Monte Carlo step. 

The vertical dimensions of the fundamental cell are de- 
termined by the minimum number of layers which gives 
accurate energies. The height of the unit cell A'layersC^ 



must be at least several times the mean in-layer inter- 
particle spacing; otherwise the computed energy will be 
artificially high due to particles iViaycrs + 1 layers apart 
which lie directly above each other. At very small d this 
limits the accuracy of the simulations, as extremely large 
numbers of layers must be used to obtain reasonable val- 
ues. 
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